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Abstract 

First order QCD and leading QED corrections to Higgs boson pro- 
duction in the channel e~p vH^X\ iJ" bb are calculated for the 
kinematical conditions at LEP ® LHC {y/s = 1360 GeV) and the inter- 
esting mass range 80 < Mh < 150 GeV. In the DIS scheme the QCD 
corrections (not including the corrections to the branching ratio, which 
are well-known) are found to be about 1% for the total cross section and 
— 13% to —10% for the observable cross section as defined by appropriate 
cuts. The latter results depend on the definition of these cuts. The QED 
corrections amount to about —5%. 



^ This work was supported in parts by the Alexander-von-Humboldt Stiftung and 
German BMFT under contract number 055ZT91I-01 



1 Introduction 



The search for the Higgs boson, or physics which replaces the Higgs boson in case it does 
not exist as a fundamental particle, is one of the main motivations for the construction 
of supercolliders. While a light Higgs boson with a mass Mh < Mw will most likely 
be found or excluded at LEP, the mass range from Mh ~ 140 GeV up to about 1 TeV 
can be completely covered at LHC and SSC. However, the intermediate region from 
Mh ~ 80 GeV to 140 GeV is a very difficult one for searches in pp collisions . In this 
mass range, the Higgs boson decays dominantly into bb quark pairs, a channel which is 
difficult to isolate because of the huge QCD background. The rare decay 77, on 

the other hand, provides a more favourable signature, but at the cost of a substantial 
reduction in rate 1^, |2) . Although it appears feasible to detect some signal in the two 
photon channel, this task is experimentally very demanding. Clearly, the best tool 
for the discovery of an intermediate mass Higgs boson would be a linear e^e~ collider 
at a minimum cms energy of roughly twice the LEP 200 energy ||^. Unfortunately, 
at the time when LHC and SSC are expected to start operation, a linear collider in 
the required energy range will most likely not yet be available. Therefore, it is an 
important issue to clarify whether or not the ep option which exists at LHC could help 
to investigate this very interesting mass window. 

The ep option we refer to can be realized by intersecting a 50 to 100 GeV electron 
beam from LEP with the 7.7 TeV proton beam from LHC jj]. In this mode, the ex- 
pected luminosity ranges from 2 • 10'^^cm~^s"^ at the upper end of the cms energy 
range, = 1.79 TeV, to 2 • lO^^cm-^s'^ at the lower end, ^/s = 1.26 TeV. For 
definiteness, we assume \/s — 1.36 TeV and an integrated luminosity of lfb~^ corre- 
sponding to collisions of 60 GeV electrons on 7.7 TeV protons and one year of running 
time. 

The basic electroweak processes in ep collisions are neutral current (NC) and 
charged current (CC) scattering through the exchange of virtual photons, Z and W 
bosons, respectively. Whereas the couplings of the standard Higgs scalar to the fermions 
participating in NC and CC scattering are strongly suppressed by the light fermion 
masses, the couplings to the W and Z bosons are of ordinary electromagnetic strength. 
Hence, the virtual weak bosons appearing in NC and CC scattering processes can act 
as efficient sources for the standard Higgs boson. In fact, the WW fusion process 
depicted in fig. ^ is the dominant Higgs production mechanism in ep collisions |^, 
followed by ZZ fusion, which has a roughly 5 times smaller cross section. Interestingly, 
in pp collisions the analogous mechanisms become important only for very heavy Higgs 
bosons with masses Mh > 600 GeV Q] . 

As can be already expected from the above remarks, Higgs bosons are not produced 
very frequently in ep collisions, at least in comparison to the corresponding rates in pp 
collisions. In the interesting mass range from 80 to 140 GeV and for y/s = 1.36 TeV, 
the total cross section varies from about 200 to 100 fb. One can therefore not afford 
to search in rare decay channels such as iJ" — > 77, the mode considered for hadron 
colliders, but one must try to detect the Higgs signal in the main decay channel bb 
(or for Mh > 130 GeV through H° WW*, where one of the W bosons is off-shell, 
this possibility is not pursued further here) . The main background is expected to come 
from NC production of jets (including photoproduction), multi-jet CC scattering, and 
single W, Z and quark production with subsequent hadronic decays. Although this 
background is large, it is not as overwhelming as the QCD jet background to bb 
in pp collisions. In addition, in ep collisions the H*^ —> bb decay provides several 
characteristic signatures which allow for an efficient discrimination of signal against 
background. For the WW fusion channel ep vH^q + X ^ vbbq + X, the selection 
of signal events and the suppression of the backgrounds have been investigated in 
great detail in ref . Q] . It has been demonstrated that one should be able to observe an 
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intermediate mass Higgs boson, provided flavour identification capabilities are available 
with efficiencies similar to the expected capabilities of the DELPHI detector at LEP 

This encouraging result and the importance of the physics issue call for further 
studies in order to corroborate the above conclusion. Obviously, the feasibility of flavour 
identification has to be examined in a dedicated ep detector study. Furthermore, one 
has not yet exploited all possibilities to optimize the event selection and background 
suppression. Finally, the parton level treatment in tree approximation performed in 
ref. still suffers from uncertainties due to the neglect of hadronization effects and 
higher order corrections. In particular, K-factors of 2 would be quite crucial since the 
estimated number of observable signal events and the signal to background ratio do 
not leave much room for losses. Therefore, it is necessary to investigate the effects 
of higher order corrections to signal and background. One cannot rely on qualitative 
arguments and guesses based on existing calculations for production processes of other 
heavy particles ||9|] since the Higgs signal considered here is defined by a set of nontrivial 
cuts on several kinematic variables. In this paper, we calculate the 0{as) QCD and 
leading QED corrections to the WW fusion process in the zero- width approximation for 
the Higgs boson, both for the fully integrated cross section as well as for the observable 
cross section in the channel — ^ hh after application of the kinematical cuts suggested 
in ref. . 

The paper is organized as follows. In section || we recalculate the tree level cross 
section for Higgs production via WW fusion and examine the variation with different 
sets of structure functions. We also show interesting differential distributions for the 
channel — > bb and briefly review the selection cuts used in r ef. 0. Section | 
is devoted to the QCD corrections. Results are presented for the fully integrated 
production cross section, including the QCD corrected branching ratio for iJ" — s- hb. 
The main point, however, is the calculation of the QCD corrections to the cross section 
times branching ratio in the three jet channel in the presence of the kinematical cuts of 
ref. 01 . The effects of initial state QED bremsstrahlung are then obtained in section ^ 
using the structure function method. In section H we summarize our conclusions. 



2 The Born Approximation 

The lowest order production cross section for e^p uH'^X is obtained from the 
diagram shown in fig. ^ The differential cross section is given by|^ 

i 

1 g^M'^ ^ r 

= ^ {Ql + M^^mi + AP^Y g [{PwMp'.)m{-. Ql) 

+ {pip'2){p'iP2)Ux, Ql)] dPS^[eq ^ vH\) dx (1) 



where the matrix element squared includes an average of the spins of incoming particles 
and a sum over the spins of outgoing particles. In the above formula pi{p'i) and ^2(^2) 
are the incoming (outgoing) lepton and quark momenta, respectively, x is the fraction 
of the proton momentum carried by the incoming partons, s the cms energy squared 
and — —{pi — p'j)"^ . Ui{x,Q^) are the distribution functions of u— and c— quarks, 
di{x,Q^) of d— and s— antiquarks. The contributions from bottom and top quark 
densities can be neglected. Also we have put off-diagonal elements of the CKM-niatrix 



^For e+p-scattering the distributions Ui, di have to be replaced by di. 
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involving the top quark to zero. The SU(2) gauge couphng g is given hy g = e/ sinOw, 
e being the electromagnetic coupling constant and 9w being the weak mixing angle. 
Furthermore, we apply the zero-width approximation, which is a valid simplification 
for Mh < 150 GeV where Th/Mh ^ 10"* 0, [l|- The phase space element dPSn is 
given by 

dPS^ = (27r)(*"3") S^'Hpin - Ep.) n f^. (2) 
The physical boundaries in x are Mfjjs < a: < 1. The integrated cross section afj^ was 



checked against earlier calculations ||, ^. Throughout this paper we have taken the 
following values for the electroweak parameters: Mw = 80.6 GeV, sin^ 9w = 0.230, and 
e^/47r = 1/128.5. The quark distributions in eq. |l| are parametrized as dependent 
functions corresponding to the solutions of the Altarelli-Parisi equations. 

In fig. H we show ctj^] and ctj^] x Br{H'^ — > 66). The branching ratio is taken from 
ref. ||l^ and includes integrated QCD corrections up to 0{a'^). The fast decrease of 
the cross section times branching ratio in the bb channel at higher Higgs boson masses 
is due to the decay WW* which becomes dominant as Mh exceeds 120 GeV. 

Concerning the dependence of the cross section on the parametrization of the quark 
distribution functions the following remarks may suffice. We use the parametrization 
by Morfin and Tung [|2| set 2 assuming the MS or DIS scheme ||l^ for the definition of 
the parton densities. The resulting CTj^j for the two schemes differ by less than 3%, as 
can be seen from fig. ^ We have also compared a^^] for other parton parametrizations. 
The deviations for HMRS (MS, set 2) are also smaller than 3%. Larger differences 
of 0(10%) are obtained for older parametrizations such as DFLM (DIS) |l5j and DOl 
(LO) 111. 

Eq. Il| was parametrized in terms of the individual parton densities directly. How- 
ever, for the understanding of the size of the QCD corrections to be discussed in the 
next section, it is interesting to know which of the usual deep inelastic structure func- 
tions, that is 

2 

2xFi°\x, Ql) = F^°\x, Ql) = 2xY,[Mx, QD + d^{x, QD] (3) 

i=l 

and 

2 

xF^°\x,Ql) = 2xY}n^{x,Ql) ~ d,{x,Ql)], (4) 
1=1 

yields the dominant part of the cross section. For this purpose, eq. |l| may be rewritten 
as 



with 



dPS^ dx 
(5) 



= \[{j>IP2){AP2) ± (P1P2)(P'1P2)]. (6) 

F!f\x,Q'^) contributes about 98% to 

'^tot- This has two reasons: b ^ 6+ in a wide 
range of the phase space|^ and xF^'^ < F^^ . 

To illustrate some details of the VFH^— fusion process differential distributions for 
X — — g|/2(Pg2) [P being the proton momentum), Q\ = — and the transverse 
momentum p± and rapidity rj of the Higgs boson arc shown in fig. ^. The figures refer 

^For ph 0, this corresponds to deep inelastic scattering at small y. 
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to a Higgs boson mass of 100 GeV. As one can see, the production involves large x 
and Q2, typically x > 0.01 and Q2 ^ 100 GeV^ . We thus do not expect screening 
corrections to the parton distributions to influence the cross section significantly]^ 
The background to the vbbX final state was investigated in ref. in detail. It re- 
ceives contributions from four types of reactions: photo(and NC)production of multijet 
events, CC multijet production, single W and Z production with subsequent hadronic 
decay, and top production. This background was shown to be manageable after apply- 
ing the following selection criteria: 

• the final state should contain three jets, defined here as partons with transverse 
momentum p± > 20 GeV, rapidity \ri\ < 4.5 and a distance AR > 1 in the r]-(j) 
plane; 

• the event should have a large missing transverse momentum, p™''"* > 20 GeV, 
and a large total transverse energy, E± > 100 GeV; 

• the Higgs signal should be contained in the invariant mass distribution of the 
two jets lowest in rapidity. 

These cuts leave from the total cross section in the bb channel, a[°lxBr{H" -> bb), the 
fraction denoted by (7^°^ in figure |[ In this sense, ct^°|. can be considered the observable 
cross section in this channel. If one furthermore assumes good flavour identification 
and a 10 GeV mass resolution the background is reduced to a level which would enable 
detection of the Higgs boson between 80 and 140 GeV 0. 



3 QCD Corrections 

The QCD corrections to e^p vH^X\H^ bb consist of three contributions: the 
corrections to the W-quark vertex W*q q' , corrections to the decay vertex iJ° bb, 
and gluon exchange between the final state 6-quarks and the interacting parton at the 
hadronic side. To 0{as) the latter corrections vanish because of the colour structure. 

We first consider the total integrated corrections. They only need to be calculated 
for the production vertex, as the integrated corrections to the decay are known jl^, [l^ 
and have already been taken into account in the branching ratio. Then we discuss the 
QCD corrections to the observable cross section, which strongly depend on the cuts 
applied to isolate the signal. In this case there is also a negative correction to the decay 
width. 



3.1 Corrections to the total cross section 

The 0{as) QCD corrections to the total production cross section are described by the 
diagrams of fig. ^. Generally, the differential cross section for ep —i- vH'^X can be 
written as 

^ ^s iQuulnZi^M^^Y + + ^^^^ 

where pc and pv denote the parity conserving and violating parts of the Icptonic [L^^'^) 
and hadronic (W^^) tensors, respectively. The hadronic tensors IV^j, are the same 
as those which appear in deep inelastic charged current scattering with = — 

^We verified this using the KMRS distributions (MS, B-5) pi. 
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(cf . ^ |T| ) . For the individual contributions to cq. ^ one obtains^ 

LM^-WPc.L ^ h^TL{x,Ql)/x (8) 

i|^:W^P' = h+T^{x,Ql)/x (9) 

- h-T:i{x,Ql). (10) 

Formally, &^ can be taken from eq. |^ with p2 ~ xP, p'2 ^ P2 + Q2 not necessarily 
corresponding to particle momenta (only at lowest order). The coefBcient is given 

by 

b^{pi,p[,P2,p'2) = - {(Pig2)(p'ig2) - iPlp'i){p2q2)} ■ (11) 

Here, the structure functions arc defined by 

J^..{x,Q') = (x,Q^) + ^^F,('\x,Q') + 0{al) (12) 

with fP{x,Q^) = and F^"^{x,Q'^) as given in eqs (|) and (|). For a^(g2) the 
leading order expression for five active flavours is used. The value of the QCD scale A 
is chosen in accordance with the parton distribution functions.^ 

For the regularisation of the coUinear divergences a factorisation scheme has to be 
chosen. Because we consider a deep inelastic process the choice of the DIS scheme 
jl^ appears to be most natural.^ Different corrections would be obtained using other 
schemes, e.g. the MS scheme. 

In the DIS scheme the structure function 

J^2{x,Q')^F^°\x,Q') (13) 

is preserved to all orders and corrections occur to the terms L^^c^Jlu'^ ip^W^"'^ 
only. In lowest order, eqs |7|-[l0| are seen to lead to da''^^ as given in eq. |. The QCD 
diagrams of fig. ^, integrated over all phase space, yield 



{x,Ql)/x + b-F^'\x,Ql)] (14) 



where 



F, 



\x,Q') = X r -{fl{z)-F^"\x/z,Q') + 4fS{z)G{x/z,Q')} (15) 

Fi^\x,Q^) = -Cp [' -{l + z)F^"\x/z,Q'). (16) 

Here, fl{z) = 2Cfz and /f (z) = 8Trz{1~z), with Cp = 4/3 and Tr = 1/2. G{x,Q^) 
denotes the gluon density. The scale in both the distribution functions and as is taken 
to be — —q2, in accordance with the definition of the DIS scheme. 

In the limit pn 0, eqs |8|-|ic| yieTd the familiar expressions for the deep inelastic scattering 
cross section: L^^Wp^ -{sx/4:fy^TL{x,Q^)/x, Lf^^Wp^ {sx/4:fY+T2ix,Q^)/x and 
L';;^WP''' ^ isx/-ifY.T3{x,Q^), with g2 = Q? = Qi and y± = 1 ± (1 - yf. 

^ Since we are calculating the finite terms of the order Os corrections, next to leading order 
parton densities and as{Q^) should be used. However, as the final corrections turn out to be 
small the difference is negligible. 

®Note that for higher order corrections one can choose the schemes to renormalize the /3 
function and the coUinear divergences separately. 
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The QCD corrections to the total cross section obtained from eq. |l^ are displayed 
in fig. ||. We show separately the quark and gluon contributions from cr£ and 

a^, and the contribution cr| from F^^\ In the mass range 50 < Mh < 150 GeV the 
correction is dominated by the term, to which quarks and gluons contribute about 
equally. The cTg^'' correction is negative and amounts to less than 10% of aj^. The total 
correction alll{QCD) adds up to less than 1% of the Born cross section. Using for 
the parametrization of the parton densities DOl ||l6| instead of MT ||lj] the relative 
correction 5^^^{QCD) = alll{QCD)/al'^l does not change significantly. This reflects 
the similar shape of the parton densities in the region in {x,Q'^) which contributes most 
to the cross section. 

As a consequence of the above result we remark that the integral 0{as) correction 
to Higgs production in pp collisions proceeding via WW fusion is of 0(2%). 



3.2 Inclusion of kinematical cuts 

3.2.1 Corrections to the production cross section 

A Higgs signal in the bb channel can only be observed after applying suitable cuts (see 
section 1^). These cuts strongly influence the QCD corrections, as exactly three observ- 
able jets are demanded. This means that there is an upper bound on the angle and 
energy of the extra gluon or quark with respect to the other three jets. These bounds 
can only be implemented numerically. We have chosen to regulate the divergences as 
follows: the ultra violet divergences using the usual dimensional regularization, the 
infra red divergences with a small gluon mass A, and the coUinear divergences (mass 
singularities) with a small quark mass m 3> A. This way all phase space integrals can 
be evaluated numerically in four dimensions. This mass regularisation scheme has been 
described in ref. [2^ . 

The total corrections can then be written as a sum of four 4-dimensional integrals: 

(QCD) = + + a« ^ + a« „*e.- (17) 

The first two originate in virtual and real gluon radiation, cr^l^g is the contribution 
from initial state gluons, while the counter terms arise from the renormalization of the 
distribution functions and are needed to cancel the coUinear divergences. 
The virtual and soft contribution can be written in the form 



, 2 

i—1 

[iPlP2){p'lP'2Wi^Lt+soft i^^ Ql) 



2xs {Ql+M^Y{Ql + M^Yy 



HPiP'2){p'iP2)d^^'Lt+softi^^Ql)] dPSsdx (18) 



i,virt-\-soft 

virt+soft 



where the functions q^j^lij-t+soft ^'^^ given in the appendix, eq. (A. 7). They still depend 



on the quark mass to and the soft cutoff A = (pJ, + fc)mm' where k is the gluon 
momentum. The expressions for da'\^^^^ and da^l^g are also provided in the appendix 



by eqs ( A.10| ) and ( A. 11 ), respectively. 



Because we decided to calculate the QCD corrections in the DIS scheme, demanding 
eq. ( |l3| ) generates a counter term to the quark distribution functions: 

(19) 
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The various contributions to ( p^ ) were derived in several previous calculations ]l9|, pO| . 
I2III using different renormalization methods. They are listed in the appendix (eqs ( A.1)- 



(A.4)) for the regularization method used in the present calculation. The ultra violet 
and infra red singularities as well as the dependence on the parameter A separating 
the soft and hard part of the gluon Bremsstrahlung terms cancel in the sum eq. ( p9| ) 
and only the mass singularity, that is the logarithmic dependence on m, remains. One 
obtains 



^i, counter 



2 



(20) 



where the functions Pqq and Pqc are given in eqs (A_^) and (A^) respectively. Using 
eq. ( po| ) one finally has 



da, 



(1) 



-i 

2xs ^ 

i—1 



\M 



(0)|2 (1) 

I I ^i, counter 



{x,Ql)dPS3dx 



2x3 iQl + M^)^{Q^, + M^) 



2 

2E 

4 = 1 



[iPlP2)ip'lP'2)u'i%unteri^^Q2) 



+ {PlP2)ip'lP2)d[ 



i, counter 



dPSs dx. 



(21) 



In the practical calculation we combine the counter term given in eq. (|l]) with the 
virtual and soft contribution, eq. (p^. Using eq. ( pj| ) one sees that the virtual and 
soft parts cancel, and only the hard and glue part of the counter terms needs to be 
evaluated. 

It turns out that the computational procedure outlined in ref. |^ (straight eval- 
uation of the individual integrals using Monte Carlo techniques) is not suitable here, 
since the corrections are very small. The problem is that we still have large logarithms 
that are subtracted after the numerical integration over the phase space. The unavoid- 
able inaccuracy in this integration (usually 0(10^'^)) then gives a very large error in 
the result. This problem was avoided by performing the subtraction of the coUinear 



singularity in the hard Bremsstrahlung and glue contributions, eqs A. 10 and A. 11 



respectively, and in the counter terms before integrating over the rest of phase space. 
One starts by rewriting the order of these integrals:^ 



hard 



counter jhard 



dz- 

dz 

2zs 

2 



(dPS^{eq ^ vH\g) e{s - A) ^ \M^'larf<l^{^, QI) 

i=i 

Jds0{s-A) jdPS3{eq^vH\*)^ JdPS2{q* ^ qg) 



X 



E 

i=l 



as 
'2tt 



dx 



1 

2X3 



1-5 



d^ 



xqi 



7 (1) 

This is for the quark contribution. The procedure for cr^^^ is analogous, 



glue 



(22) 
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- / rf-^ / T ^(1 - ^ - S) ldPS,{eq ^ .H\) Y: \M[y U, %) 

T^Jt) ^ZSJq t; J .^^ \ 1^2' 

xq,{z,Ql) (23) 



with s = {p2 + kY and (5 = (A — m?')/Q\. Next, we also introduce in the counter term 
a gluon momentum k — P2{^ ^ C)/C ^-^ foUows: 

jdPS^{eq{p2)^vH\)x Jd^k5('\k^^p2)e{l-^-d) 



e 



dPSiieqips) ^ i.i7"gg)(2^)32fco -J^'^fc - (1 - Ops) 0(1 ~ ^ ~ S) 



Jds JdPS3{eq{p3) ^ i^H^*) X J ^ j dPS^iq* qg) (2^)22fco^^'Hfc-(l-OP3) 
X 61(1 -^-(5) 



Ids JdPS^ieqips) ^ iyH\*) ^^Oi-^ ~ ^Y^^) (24) 



where p2 — xP and P3 — P2 + k ^ zP, and terms of ©(m^/Qj) are neglected. The 
2 — !■ 3 phase space integrals in the hard Bremsstrahlung and counter term now have 
the same kinematical boundaries. This means that the subtraction can be performed 
for each point in this phase space. The integrals over s are also seen to be identical in 
the limit A — > 0. What remains are integrals over cos 6* and for the Bremsstrahlung 
term (cos0, are the angle variables in the q* cms), which do not have a counterpart 
in the counter term. 



In the integration of the hard and glue contributions, eqs A. 10 and A. 11, it will be 
necessary to treat the single and double pole terms separately. The double pole terms 
are non-divergent (the numerator is proportional to m^), but give a finite contribution. 
They have a malicious cos 9 dependence, but can be integrated numerically with a 
suitable mapping. 

In the rest, the integral over cos 9 gives an exact cancellation of the single pole terms 
in the coUinear limit cos 9 —^ 1 when the radiative te rms are sy mme trized in cos 9, and 



the collinear logarithm in the counter term (see eqs ( A. 5 ) and ( A. 6 )) is represented by 



ln(^(l - 0''mVQl) = / :^»^2.o2^-1 (25) 



d cos 9 

cos0-(i-e(i-O"™VQl)-' 

with n = ±1 for the quark and gluon contributions, respectively. Now the sum of 
the radiative contribution, excluding the double pole terms, and the counterterms is 
finite for all values of cos 9. As the boundaries on all integrals coincide this sum can 
be integrated numerically. 

For clarity we repeat that as the virtual and soft contributions were seen to cancel 
against their respective counter terms, only the above integrals, which combine the 
hard Bremsstrahlung and gluonic contribution with their respective counter terms, 
have to be evaluated numerically. The independence of the result on the soft cut-off A 
and the quark mass m was verified by varying the value used in the computation. The 
numerical accuracy now suffices: without cuts we find agreement with the analytical 
calculation presented in the previous section within the statistical accuracy (1-3%). 

When applying the cuts one has to define the way to combine the four jets produced 
in eg —> vH^qg and eg —> vH^qq with — s- hh into three jets. The algorithm used 
was to combine the two partons with minimal distance in the rj-cj) plane if this distance 
was less than 1, and else to require that one of the partons did not pass the (p_l,t?) 
cuts. The resulting 3-jet event was then subjected to the cuts described in section H. 
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In the coUinear and soft limits this is seen to be equivalent to the cuts applied to the 
Born approximation.^ 

With the experimental cuts applied a numerical accuracy of about 10% of the 
corrections is achieved. The reduction in accuracy is caused by small regions in phase 
space, close to the collinear limit cos = 1, where the three jet signal is cut but the four 
jet survives, or vice versa. Of course at the collinear limit the kinematical cuts are the 
same, but for any finite angle between the quark and gluon they differ. The resulting 
wedge-shaped regions have the original logarithmic singularity, and hence give a finite 
but small contribution. Unfortunately present-day numerical integration methods 
do not handle such highly discontinuous functions in a high-dimensional phase space 
very well. 

The results of this calculation are summarized in fig. ^. One can see that the effect of 
the cuts is to make the corrections negative and larger. This is to be expected, as the re- 
quirement to observe exactly three jets cuts away part of the (positive) Bremsstrahlung 
terms. The almost exact cancellation which keeps the integrated corrections so small is 
now violated, thus are the relative corrections larger with cuts taken into consideration. 
The gluonic contribution is much smaller as the gluonic component tends to be softer, 
which reduces the probability that four jets are observed. 



3.2.2 Corrections to the decay 

The 0{as) correction F*^^) to the decay width of the Higgs boson into 6-quarks has long 
been known . The diagrams for this correction are shown in fig. |^. Recently also the 
O(a^) contributions F^^^ were evaluated |2|]. In the mass range 80 < Mh < 150 GeV 
and adopting the MS scheme, the relative corrections F^^+^^/F^^^ vary from —28% to 
—43%. In the applications considered so far in this paper only the branching ratio 
Br{H'^ bb) enters. This is affected less than the bb partial width, at least as long 
as H'^ bb is the dominant decay mode. Numerically, we derive from ref. ||ll| that 
-3.4% > Br(i+2) /bA"'> > -36% in the mass range considered.| 

However, these corrections also change in the presence of kinematical cuts. For 
the moment, let us consider the 0{as) corrections only. A part of these corrections 
involves the emission of a hard gluon, 7?" — > bbg, which can change the configuration of 
the bottom jets considerably. These events should not all be included in the observable 
cross section as we demand exactly three jets, of which two originate in the bottom 
quarks. The difference must be considered an 0{as) correction to the results given 
in section ^ which already include the QCD corrected branching ratio Br{H^ —f bb). 
This correction is scheme independent as only the hard radiation part is involved. It 
is also expected to be relatively small, as the main effect of the integrated corrections 
is to introduce a running b quark mass in the Yukawa coupling. These terms are 
of order ln(M^/m^); they arise from the renormalisation of the bbH coupling and 
do not contribute here. Furthermore, the effect of the cuts clearly makes the 0{as) 
contribution negative. 

The 0{as) correction to cr^°^ described above is given by the difference 

^decay " ^fcbg ^bb y'^^) 

where 



2 

"'^^ = jdx jdPS^{eq^vH\)Y,\Mf^\,{x,Ql) 



*The same algorithm was used in the calculation presented in [0] to estimate the contribu- 
tion to the background from four jet events. 

^The QED and electroweak corrections to r(_ff'' — » bb) were found in ref. to vary 
between —0.6% and —1.7% for 50 < Mh < 150 GeV. They have not been included here. 
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^ JdPSsiH' ^ bbg)\MHO^bbg\'0{4-iet cuts) / Tt^t (27) 



is the contribution from hard gluon radiation subject to the appropriate 4-jet cuts, 
while 

- U Lpc„..„^.,^o„,^,,.(o),.„.,.^.Jrf^5.(i/"^6&)^?(3-jetcuts) 



bb 



Jdx jdPS,{eq^.H\)Y:Mf\\.Ml) ^dPS,{H^-.bb) 
2M^ JdPSsiH'' bbg)\M 

H°^bbg I' / Tt,t (28) 

is the corresponding integrated contribution corrected for the effect of the 3-jet cuts. 
Ttot is the total width of the Higgs boson. Obviously in the absence of cuts (^^^^^^y ~ ^■ 
The decay matrix element, summed over all spins, \MH°^bbg\'^i does not depend 
on the orientation of the b quark in the Higgs boson cms system; only the 4-jet cuts 
will be influenced by it. We thus extract the integral over these angles J dVlb from the 
3-body decay phase space integral. The 2-body phase space integral contains a similar 
angle, which enters only in the 3-jet cuts. Identifying these two angles we obtain 



decay 



2 

(dx fdPSsieq ^ z^ffOg)^ |Alf)|29,(x,Q2) 

i=l 

" 2kji¥l jd^dE,dE,\M,.^,,X /Vtot 

X jd^lb {6l(4-jet cuts) - 6l(3-jet cuts)} . (29) 

where (f) is the azimuthal angle of the b momentum with respect to the b momentum. In 
a Monte Carlo integration this means that for each Bremsstrahlung event one constructs 
a decay into only bb with the b quark momentum in the same direction (in the Higgs 
boson cms). One examines whether this alternative kinematical configuration passes 
the 3-jet cuts: the non-zero region is the one in which the 3-jet configuration passes its 
cuts, but the 4-jet does not, or vice versa. This will not happen in the limit of soft or 
coUinear gluon emission, so the result is free of singularities and ln(M^/m^) terms. 

Finally, we can include important 0{a'j.) effects by reintroducing the 2-loop inte- 
grated branching ratio as 

\MHO^bba? _ \MHO^bbg? „ (0+1+2) .on^ 

f;;^ - r(o+i+2) ■ ^"^"^ 

The fraction on the r.h.s. of eq. ^ is now almost independent of the (running) b quark 
mass, the main effect of which has been absorbed in the branching ratio. Effectively, 
we thus use the two loop expression for the fefeiJ" coupling (which is proportional to 
the b quark mass) everywhere, while considering the effects of kinematical cuts only up 
toO(a,). 

The expression for |A^//o^f,5gP used agrees with previous calculations It 



is given in the appendix (eq. (A. 12)). The resulting correction cr^g|,Q„ is shown in 



- -, , . decay 

fig. o. The numerical result is not sensitive to the precise value of nn, in the range 
3 < mf, < 8 GeV. The correction cr^^l,^^ to the decay is larger than the correction cTp^^^ 
to the production process as the probability to detect a fourth jet radiated from a final 
state 6-quark is higher. The rise for low Higgs mass results from the high probability 
that the b and b jet become indistinguishable when a hard gluon is emitted in a boosted 
system. 
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It should be noted that the size of these corrections critically depends on the way 
in which the jets are reconstructed. Here we employed a rather simple method, which 
could be improved upon by using more sophisticated algorithms. The Monte Carlo 
approach given here allows such studies to be performed in a straightforward way. 



4 QED Corrections 

The dominant contributions to the QED corrections for ep scattering at high energies 
can be described in the leading logarithmic approximation ||2^, ^ . For charged current 
processes, only initial state bremsstrahlung (see fig. ^) from the incoming electron line 
has to be considered, since radiation from the quark line merely leads to a small extra 
contribution of 0{3ae^/4:as) to the QCD evolution. In this approximation, the QED 
correction to the Higgs production cross section is given by 

a«(Qi?i?)^£|dzi±^ln(£^){a(")(zpi)--^°Hpi)} (31) 

where a = 1/137. The factorization scale /i is related to the fc^— integral of the radiated 
photon, ^"^ Q\Ki s. We have chosen jj.^ = s (sec ref. |p7|). 

In fig. ^ the QED correction a^^^{QED) to the total production cross section is 
shown as a function of Mh using the parametrization MT (DIS, set 2) |lj]. In the 
mass range considered, it amounts to —4% to —5% for the choice /i^ = s. The QED 
radiative corrections are rather small because of the small ratio Mfj/s < 0.012. A 
similar behaviour was recently observed for the correction to heavy flavour production 
in neutral current ep— scattering where for small values of Auiq/s only small neg- 
ative corrections are obtained. The correction to the observable cross section after cuts 
is also shown in fig. |. The relative correction (T^obsiQ^D) / '^^ol = -4.5% to -5.8%.[3 
These results are only weakly influenced by the choice of the parton distributions. 



5 Conclusions 

We calculated the 0{as) QCD corrections to Higgs boson production in e^p collisions 
via the WW fusion process followed by the decay bb in the zero-width approx- 

imation. Results are obtained for the fully integrated cross section, as well as for the 
observable cross section after application of the kinematical cuts suggested in ref. 0] . 
We find surprisingly small corrections of the order of a percent and less for the uncut- 
ted cross section, and larger and negative corrections when the cuts are included. Our 
results are at variance with the expectation expressed in ref. Q. 

We conclude that the QCD corrections that have to be added, once the known 
corrections to the branching ratio are included, in order to obtain the complete 0{as) 
corrections to e^p vH'^X\H^ — > bb, are in fact very small; the only sizeable effects 
originate in the selection cuts which mainly affect the hard gluon radiation. The 
following factors contribute to make the corrections this small: 

• The cross section is dominated by F2, which is not renormalized in the DIS 
scheme. The difference between the DIS and MS structure functions suggest 
that the corrections would change by a few percent when changing to the MS 
scheme. 

• The scale is large everywhere [Q^ > 100 GeV^). This means that as is small 
{as/2TT « 0.02). 

^"The leading QED corrections to the W and top backgrounds were also calculated and 
found to be comparable to the corrections to the signal. 
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• No large logarithms enhance the corrections, as these have already been absorbed 
into the evolution of the distribution functions. 

The main effect due to the cuts occurs at the decay vertex, where the extra gluon is 
either detected as a fourth jet or causes the b jets to be too close together. These effects 
however depend on the algorithms used to define the jets. 

In addition, we also calculated the effect due to initial state QED bremsstrahlung 
from the electron. This correction decreases the observable cross section by roughly 5%. 
The relative QCD and QED corrections to the observable cross section are summarized 
in fig. One can conclude from our results that the analysis performed in ref. 
is not changed significantly by higher order effects as far as the signal is concerned. 
What remains to be done is a corresponding calculation for the background processes. 

Acknowledgement. We would hke to thank J. A. M. Vermaseren, K. J. F. Gae- 
mers and D. Zeppenfeld for discussions. 
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A Appendix 



In this appendix we list explicit expressions used in the calculation of the 0{as) QCD 



corrections to the observable cross section as discussed in section 3.2 



In the regularization scheme adopted in this calculation, the 0{as) modifications 
of the quark distributions due to the virtual, soft and hard Bremsstrahlung terms and 
the gluon contribution are given by 

-l2 



a^^^ (x 



3_ 2^ 
2^ 



In 



2 In 



A2 



1 + ln 



(A.l) 
A 



27r 



l-<5 



Pn. 



with A = (p^ + (5 = (A - m'^)/Q'^ and 



i+e 



In 



= -<;(i-2e + 2ain 



1-e 



-8^ + 8^-1 



The dependence on A cancels in the sum of the soft and virtual parts: 



27r 



7 . 



A 



(A.5) 
(A.6) 



ln(A 



xq,{x,Q'^). 



Furthermore, ql^hardi^^ Q^) ^le rewritten as 



as 
2tt 



i-s 



/ Tfi \ r -L 

nl-<5 



(A, 



In the first integral the limit (5 — > can be taken because terms of 0(S) may be 
neglected. The second integral can be calculated analytically yielding 



^^F<i ^ + ^ + Iln(5 + ln"<5 + 21n( ^ 



ln(5+ - 
4 



qi{x,Q^). 



(A.9) 



In the limit m ^ 0, the contribution ( [A.9| ) just cancels the term ( |A.7| ) in the sum 

qi^virt-^rsoft ~t~ qi^hard- 

The differential cross sections for hard gluon Bremsstrahlung for an incoming quark 
or antiquark and the gluon contribution were derived using FORM [^0[ and are given 

by 



d^i'L = ^^T.\^^iard\'<l^i^,Ql)0iS-^)dPS4 



dz 



} 



m 
(A.2) 

(A.3) 
(A.4) 



1]} 

(A.7) 
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^ AirasCF^l [~{kpi){p[p'2){p3P2)/{kp3)/{kp2. 



2zs {Ql + M^y{Ql+M^r 



+ (fcpi)(p'iPi)/(fcp3) + ikp\)ipiP3)iP3P2)/ikp3)/ikp'2) + ikp\){piP3)/ikp'2) 

- iPlP3)ip'lP3)/ikp3) + 2ipiP3)ip'iP2)iP3P2)/ikp3)/ikp2) ~ iPlP3){p'iP2) / ikp3) 
+ {PlP3WlP'2)/{kp'2) + {PlP'2){p'lP'2)l{kp'2) + mHkp^){p',p'^)/{kp3f 

- m\kp',){j>^p3)/{kp'^f ~ m\p,p3){p',p'2)/{kp3? - m'{piP3WiP'2)l{kp'2f] 

X u,(z, Ql) + [pi <-> p\] d^{z, Q2)| 0{s - A) dPS^eq vH^qg) dz, (A.IO) 

dc^iL = ^JMllfG{z,Ql)dPS,dz 

,4^,-2 +(fcpi)(Vi)/(fcp3) 



2zs(g2+M2,)2(Q2+^2^)2 «2 

+ 2{kp,){kp'Mp'^)/{kp3)/{p3p'2) - {kpi){kp'2}{p[p3)/{kp3)/{P3P2) 

- ikpi){p[p'2)/{kp3) - ikpi){p[p'2)/iP3P2) + {kpiKp[p3)/{p3P2) 

- {kp'2){piP3){PlP2) / {kp3) / {P3P2) + {PlP3){PlP2)/{kp3) + {PlP'2) {PlP'2) / {P3P2) 

+ m'{kpMp3)l{p3p'2f - m\kpMp'2)/{kp3f - m\kpMp'2)/{P3p'2? 
+ m\p^P3){p[p'2)/{kp3?'^G{z,Ql)dPS^{eg^uH\q)dz. (A.ll) 

Here k denotes the momentum of the extra outgoing parton and p3 = zP is the 
momentum of the incoming parton. The decay —> bb is again inserted when the 
observable cross section is calculated. 

The matrix element squared for the hard Bremsstrahlung part of the decay width 



used in section 3.2.2 is given by 

\MHO^bbg? - 4^asCF|+8 + 4(pbfc)/(pfcfc) + 8{pbPi) / (pbk) + KPbPl) I {Plk) 

+ HPbPbf/{Pbk)/{Pbk) + 4(pbk)/{pbk) - 4ml{pbk)/{pbkf - AmlipbPi) / {pbkf 

- 'iml{pbPb)/{Pbkf - 8ml{pbp-b)/{pbk)/{Pbk) - 4ml/ (pbk) - '^'^1 / {pbkf (pik) 

- 4ml/{pbk) + 4mt/{pbkf + Aml/{pbkfy (A.12) 
The two- loop corrected decay width for bb is derived in ref . : 

p(o+i+2) ^ _^ b y H J ^ ^ 303^^ ^ 2.953a2), (A.13) 

3z7r Mjj 

with the running b quark mass rhb and coupling constant as taken at the scale Mfj. 
As stated in the text, the magnitude of the 0{as) corrections discussed in this paper 
does not depend significantly on the value chosen for this mass. 
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B Figures 

Figure 1: The Born diagram for Higgs boson production via WW fusion. 

Figure 2: Born cross sections using the parton distributions by Morfin and Tung 
(MT2). Long-dashed line: 4'^] for MT2(MS); full line: (t|°] for MT2(DIS); 
dash-dotted line: 4°] x Br{H^ bb) for MT2(DIS); short-dashed line: cr^J^ 
as defined in the text. 

Figure 3: Differential distributions in x, Q^, the transverse momentum p±^ and 
rapidity r] of the Higgs boson for Mh = 100 GeV using the Born approximation. 
The dashed lines give the distributions after application of the selection cuts 
proposed in ref. [^. 
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Figure 4: The 0{as) diagrams contributing to the QCD corrections at the 
hadronic vertex. 



Figure 5: The 0{as) QCD corrections to atot in the DIS scheme. Full line: 
alll{QCD), long-dashed line: o"^, dash-dotted line: a^, and short-dashed line: 



4- 



Figure 6: The 0{as) QCD corrections including the cuts explained in section |2|. 
The corrections to the production process ep — > vH^X from quarks and glu- 
ons and to the decay process bb are shown separately. The total sum 

a^Ji^^{QC D) is the appropriate 0{as) correction to the observable cross section 

'^obs (which already includes the QCD corrected Br{H^ — > bb)). The error bars 
indicate the precision of the Monte Carlo calculation. 



Figure 7: The 0{as) diagrams contributing to the QCD corrections to the decay 
vertex. 



Figure 8: Diagram representing the leading QED correction to charged current 
ep processes. 



Figure 9: The leading QED corrections to the total (dashed line) and observable 
(full line) cross sections for e~p — > uH^X; bb (using fi'^ = s). 



Figure 10: The relative corrections to the observable cross section 

for e~p — > uH^X; — > bb, including the cuts given in section ^. The QCD 
and the leading QED corrections are shown separately. The error bars indicate 
the statistical precision of the Monte Carlo calculation. 
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